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1 Introduction 

This study gives the development of a new technique for analyzing images of Cryogenic 
sprays, to estimate the drop-size distribution. It has a sound mathematical basis, in 
the form of Mathematical Morphology, and we have tried to build up a formulation for 
a granulometry, starting from the elementary operations of Dilation and Erosion. An 
axiomatic foundation for granulometry has also been discussed. 

We have taken an actual LN 2 spray photograph for analysis, to illustrate the use 
of Morphological operations, culminating in a granulometry. 

2 The Malvern Mastersizer 

The Malvern mastersizer is based on the principle of laser ensemble light scattering.lt is 
currently, the most widely used instrument for estimating the droplet size-distributions 
of sprays. Ensemble scattering means that droplets cannot be measured individually; 
only a group of a large number of droplets can be sized. It is a non-imaging instrument, 
since sizing is done without forming an image of the object. 

A laser beam emerging from the source, is scattered by a collection of droplets in 
a spray, and the scattered light falls on the receiver. The receiver is an array having 
several annular detector elements. Each element senses the light scattered by droplets 
in a particular size-class. The sizes of the droplets are related to the angle of scattering. 

3 Need for an Imaging system 

In case of the Malvern Mastersizer, sufficient light must fall on the detector for it to 
sense the amount of scattering. Excessive obscuration, encountered in dense sprays 
prevents this. 

Further, in the specific case of a liquid nitrogen spray, injected into the ambient at 
supercritical temperature, vapour around the evaporating droplet leads to density gra- 
dients, which in turn affect the refractive index. This effect is called beam steering. 

As it is, Malvern Mastersizer gives an excessively large SMD (typically, 800/im 
at 450kPa).We know from the physics of droplet breakup, that this is not possi- 
ble. Moreover, studies by Chin et.al [3] indicate that the Malvern Mastersizer under- 
predicts the SMD, of dense sprays. Thus the drop size distribution given by Malvern 
Mastersizer is physically untenable. There is, therefore, a need for an imaging system, 
which can directly capture the spray on an image, which can then supply information 
about the drop sizes. 

4 Method of capturing images 

We have used LN 2 as a simulant for cryogenic propellants. The axisymmetric jet 
issues out of the plain injector, into the ambient at pressure drops between 300kPa and 
700kPa. For capturing the images a 1540-P2 strobolume pulse light generator together 
with a digital camera of resolution 300dpi is used with front-lighting. The width of the 
pulse is 50 ns at a light power output of 250 mW. The wavelength of emission is 660 
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nm with a spectral width of 27nm. Images are captured at a high film speed equivalent 
to ISO 400 on a conventional photographic film. 

5 Image Processing 

Photographs of the spray yield just a raw image, which can scarcely be of any use. 
It must be first pre-processed to a form on which Morphological operations can be 
performed. Pre-processing includes conversion to grayscale, followed by contrast en- 
hancement if required. The grayscale version is then converted to binary, using an ap- 
propriate threshold. These are standard image processing operations, and details may 
be found in [1] 

6 Mathematical Morphology 

Morphology is a term commonly used in the biological sciences for describing the form 
or structure of an organism. In image processing , by morphology, we mean the geo- 
metrical or textural features of the image. 

The techniques of morphology try to duplicate the manner in which a human being 
perceives an image. While observing a flower for example, we might first observe the 
colour, then the texture, followed by the geometry and so on, in succession not all at 
the same time. We then try to look at that aspect, which interests us the most, in 
greater detail. 

In morphological image processing(MIP) too, we select a specific characteristic of 
interest, and analyze the image with respect to that. Know ledge about the other char- 
acteristics is now irrelevant to us and, is off-loaded before we proceed. This is why in 
our experiments, the binary version has been used, though we could as well have used 
the grayscale version. 



6.1 The basic morphological operations 

Here are some basic operations which are extensively used in MIP.The book by Gon- 
zalez and Woods (Tj is an excellent introduction to Morphology. Let A be a set in 3? 2 
and x G 9ft 2 

1. Translation 

A + x = {a + x : a £ A} 

Which is to say that a point z6 A + x^3a point a G A, such that z=a+x. 

2. Minkowski addition 

For A and B C 3? 2 

A © B = UA + x, with x e B. This leads to: 

A © B = U{x + y}, such that x G A and y G B. Several properties immediately 
follow this definition: 
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Figure 1: Dilation of an image by a structuring element 



(a) A + = A. 

(b) A@x = A + xV x e$l 2 . 

(c) Commutativity, A © B = B © A. 

(d) Associativity A © (B © C) = (A © B) © C, and in particular A © [5 + x] = 
[A®B]+x 

3. Minkowski subtraction 

AQ B = HA + x, such that x <E B. With scalar multiplication being defined as 
£A = {tx : x E A} we get —A = {— x : x G A}, which is called the reflection of 
A. 

4. Erosion An alternative expression , is: 

AQ B = {x : — B + x C A}, from which the erosion is defined as: 
A © -B = S(A, B) = {x : B + x C A}. 

5. Dilation It is the same as Minkowski addition 
®(A,B) = A@B. 

Dilation and Erosion obey: 



(a) @(A,B) = 3>(B,A) 

(b) ®(A, B + x) = @{A, B) + x. 
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(c) £(A, B) = £A,B + x + x. 

(d) For a given B, and A\ C A 2 , 

i. @(A u B)c@(A 2 ,B) 

ii. C S{A 2 ,B). 



A 




B 

AqB = {x:-B+xcA} 



Figure 2: Erosion of an image 

(e) Opening Opening is a combination of dilation and erosion and is denned 
as: 

0{A,B) = @[g{A,B\B\ 

(f) Closing &(A, B) = g[9{A, -B), -B]. These operations will now be ana- 
lyzed in detail, since they were used to analyze cryogenic spray photographs 
in this exercise. 
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(-B +x) n A ^ 



Figure 3: Dilation of an image 



AE-A&S-01 



7 




Figure 4: Opening, as erosion followed by dilation 



AE-A&S-01 



8 




Figure 5: Closing operation 

6.2 Analysis of Opening and Closing 

1. <g(A,B) c = <ff(A c ,B) and 

2. 0(A, B) c = tf(A c , B) ... Duality. 

Thus, if the complement of the first operand, A is known, we can find the opening, 
given the closing and vice-versa. 

3. G{A, B) = U{B + y : B + y C A} 

which says that the opening is the union of translates of B which are contained 
within A. 

4. 'if (A, B) = n{(B + yf :(B + y)c A c } 

Closing is the intersection of the complements of certain translates of B. 

5. A point z is G 'if {A, B) (B + y) n A ^ (j) V(B + y) 

6. 0(A,B) C A 

The opening is antiextensive. 

7. A x cA 2 ^ 0(A u B) C 0(A 2 ,B) 
It is increasing. 
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8. 0[0(A,B),B] = &{A,B) 
and also idempotent. 

The closing is also increasing and idempotent, 

9. ff{A, B) D A, but is extensive. 

A set A is B-closed(or B-open) when ^(A,B) = A(or<ff(A,B) = A). 

10. A is B-closed <^ A c is B-open. 

11. This is connected with the Minkowski, addition and subtraction as: 
A is B-open<^> 3 a set E such that A = E © B and 

A is B-closed-v^ 3 a set E such that A = E B 

12. Generalized laws of idempotency for a set A which is B-open, and any other set 
F: 

G\<ff{F, B), A] = &{F, A) and 
0[0{F,A),B\ = ff{F,A) 

13. If r > s 0, with r, s, € 9ft and B is convex^ then for any A, €?(A, rB) C 6 {A, sB) 



*A set is called convex when a line joining any two of its points, lies entirely within the set. 
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7 Granulometry 

Let A be a compact set C 9ft 2 [] and t 6 SR. We, next consider \l> an image-image 
transformation(which could be the opening, on a cryogenic spray image, for instance), 
applied to A 
* : A i — ► . 

If E is a convex, generating structuring element (which is really, an image). Then {tE } is 
a family of structuring elements. The mapping t i — > 0(A,tE)\s decreasing. This means 
that^(A, tE) shall be empty * for sufficiently large t.The mapping t i — > &(A,tE)is 
known as a granulometry. Granulometric filtering is analogous to the physical operation 
of sieving. Intuitively, it is sufficient to understand that this is a kind of filtering in which 
those sections of the image which are not sufficiently large to hold the structuring 
element tE are removed from the image. 

Though this can be better appreciated, by following the algorithm, with a photo- 
graph of an actual cryogenic spray as will be seen in this sub-section, the axiomatic de- 
velopment of granulometries, given in the next sub-section, should not be skipped if one 
wishes to really understand the process. This study is on Binary images;Morphological 
operations in the gray-scale domain are discussed in [2]. 



7.1 Sieving and Granulometry on a Cryogenic Spray Photo- 
graph 

The algorithm for generating a granulometry is represented by the flowchart shown in 
Figure [6] below. It accepts the raw image A as the input and delivers a single parameter 
granulometry $ t (A) with t = 1,2, ... as the parameter. The generating structuring 
element used, E, is a unit square. 

Incrementing the value of t after each step yields a sequence of images under the 
granulometry t t— > &(A,tE). The process stops when at some stage the opening is 
empty(as it shall be; see Section [7j ir) Figure [7] shows the physical analogy for the 
operation. A sequence ^t(A) of openings is thus produced.lt can be seen that there 
will be a minimum mesh size at which nothing will be retained on the sieve. 

After a sequence of images ^t(A) , t=l,2,. . . has been generated, their complements 
w.r.t the entire image are determined. This sequence is denoted as $t(A) in Figure 



TTj which shows how the difference between successive images of this set can be used 
to get an image of droplets in a particular pixel-size class. These images are denoted 
as $ rir _i, where r denotes the higher size-class. The pixel-size classes are converted to 
length classes, knowing the resolution at which the image was captured. 

Once this is done, we just have to count the number of droplets in each such image 
Q r ^-i(A) and these values can be directly used for plotting the droplet size-distribution, 
followed by calculations for the SMD and other diameters. 



1. 
2. 
3. 



A compact set is closed and bounded. 
A closed set contains its own boundary. 

A bounded set is contained within a disk of finite radius. , centred at the origin. 
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Figure 6: Flowchart for generating a granulometry 
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(remains In sieve) 




Figure 7: Physical Analogy (Sieving) for Opening 



Figure 8: Opening by a structuring element of size 4 units ^^(A) 




Figure 9: Opening by a structuring element of size 8 units tys{A) 



Figure 10: Droplets in the pixel-size class 7-8 




4> 4 (d<4) 



*„ (4<=d<5) 



Figure 11: How the image r _i) is obtained 
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Figure [8] and Figure [9] are openings, while Figure 10 shows the droplets in the pixel-size 



class 7-8 for an image of an actual cryogenic spray. 

7.2 Axiomatic development of granulometric mappings 

Let X be the collection of Euclidean images C Then a granulometry on X is a 

family of mappings 

* t : X i — ► X, t > 0, such that 

1. %(A) C A V t > . . . ty t is antiextensive. 

2. A C B * t (A) C * t (S). . . * t is increasing. 

3. % o ^ = ^ o ^ = * max{tt{) , V t, t > 0. 

In particular, ^(^4) = AOnce we have these axioms, several results can be easily 

derived. One such result is: 

For a granulometry if r > s, ^f r (A) C ^ S {A) 

If we add the following two axioms to the existing three, we get an Euclidean 
granulometry For any t > 0: 

4. \l/t is compatible with translation. 

5. For an image A, * t = t ^(l/t A) 
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8 Conclusions 

In this report, it is seen how Mathematical Morphology can be used in processing the 
image of a spray to quantify the structures observed in it. Though this technique was 
developed, considering the need for analyzing Cryogenic sprays, it could be applied to 
storable propellants, as well. 

Morphological opening is however, a shape-distorting operation and therefore, no 
attempt must be made to relate the shapes of droplets, in the images comprising the 
granulometry, with those in the actual spray. 
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